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Abstract. We have studied the orbital ordering (OO) in Lao.5Sri.5Mn04 and its soft x-ray resonant diffrac- 
tion spectroscopic signature at the Mn L2,3 edges. We have modelled the system in second quantization 
as a small planar cluster consisting of a central Mn atom, with the first neighbouring shells of oxygen 
and Mn atoms. For the effective Hamiltonian we consider Slater-Koster parameters, charge transfer and 
electron correlation energies obtained from previous measurements on manganites. We calculate the OO 
as a function of oxygen distortion and spin correlation used as adjustable parameters. Their contribution 
as a function of temperature is clearly distinguished with a good spectroscopic agreement. 

PACS. 61.10.-i X-ray diffraction and scattering - 71.30.-|-h Metal-insulator transitions and other electronic 
transitions - 71.10.-w Theories and models of many-electron systems - 78.20.Bh Theory, models, and 
numerical simulation 



1 Introduction 



Recent soft x-ray resonant diffraction measurements on 
. half-doped manganite Lao.5Sri.5Mn04 have been inter- 
' preted in terms of orbital ordering (00) |1I2I3| . There 
, is a strong interest in understanding how 00 settles in 
■ this system as a result of the interplay between charge, 
' orbital and spin degrees of freedom. In fact, this interplay 
. can be found in many other different phenomena, such as 
' high-temperature superconductivity, colossal magnetore- 
, sistance and magnetostructural transitions 0]. 

' The Mn02 plane with its superstructure for half dop- 
ing at low temperature is shown in Figure 1 '5"6''7'8'9'10' . 
The two charge-separated Mn sites, which we shall denote 
, by Mn^+ and Mn'*"'" (although the charge separation is 
< fractional) display a checker board alternation. One can 
see ferromagnetic zig-zag chains, where the Mn*"*" sites 
form the corners and the Mn'^"'" sites are in the middle of 
the straight segments. Adjacent zig-zag chains are antifer- 
romagnetically aligned with respect to each other. There 
is a distortion of the oxygen atoms consisting of an elon- 
gation of the Mn^+-0 bonds along the zig-zag segments. 
Figure 1 shows the occupied Mn^"^ Cg orbital under the 
hypothesis of Sz^— r^/ Sx^—r"^ ordering, which is a possi- 
ble simplified way of looking at the electronic structure, 
although so far there is not yet a consensus regarding the 
correct ordering. More importantly, there is no consensus 
yet regarding the way the 00 and spin-ordering settle in 
the system as a function of temperature. 



Based on recent soft x-ray resonant diffraction experi- 
ments we will try to answer these still open questions with 
a theoretical analysis of the spectra. The soft x-ray reso- 
nant diffraction experiments that we consider have been 
done on Lao.5Sri.5Mn04 at the Mn L2,3 edges where the 
scattering factors exhibit huge variations with a strong 
dependence on the precise 3d orbital occupation and or- 
dering. 

Considering the stoichiometry in Lao.5Sri.5Mn04, the 
average occupation of the Mn 3d shell can be estimated 
at about 3.5 electrons. The soft x-ray diffraction spec- 
tra show a strong multiplet structure due to transitions 
3d" ^ 2p^3d"+^ UJj which are complicating the analy- 
sis. Previous studies [TI2..3) simplified the analysis of the 
multiplet structure by considering integer charge ordering 
(CO) with a 3d occupation of 4 electrons on the active 
sites (active in the 00 diffraction process). The scatter- 
ing factors of the active site (Mn^"*") were then calculated 
in the crystal-field approximation. Although such calcula- 
tions were able to produce some reminiscence of the ob- 
served spectra, the agreement was not good enough to 
clearly discriminate between the different hypotheses con- 
cerning the occupied orbitals. Moreover, these calculations 
were based on a local mean-field approximation and were 
inherently unable to treat the effects of spin correlation be- 
tween the different Mn sites. The approximation of integer 
occupation is contrasting with band calculations results, 
which estimate the charge separation between Mn sites to 
be about 0.3 (Ref. T?) instead of one electron. 
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Fig. 1. The Mn02 plane and its super structure in the case of 
half doping at low temperature. 

In this work we will go beyond the hmits of previ- 
ous analyses, by building a cluster around the active Mn 
site. We consider a model where the degenerate 3d shell 
of the active site is coupled to the neighbouring oxygen 
orbitals by hopping terms. On their turn, the oxygen or- 
bitals are coupled to the 3d orbitals of the inactive Mn'^+ 
sites (these for symmetry reasons inactive sites do not con- 
tribute to the 00 diffraction). For the inactive Mn sites 
we consider only one eg orbital. For this orbital we add a 
spin-dependent energy term to simulate the spin magneti- 
sation. In this model we can play with two main ingre- 
dients, i.e., the distortion of the oxygen positions (Jahn- 
Teller distortion) which rescales the hopping, and the spin 
magnetisation of the inactive sites. 

In Sec. 12 we define the Hamiltonian and in Sec. 13 we 
fix the parameters. In Sec.01we compare the experimental 
spectra to the calculated spectra obtained by choosing an 
optimal distortion. The study of the dependence of the 
calculated spectra on distortion and spin magnetisation is 
summarized in Sec. ]S\ 



2 Model Hamiltonian 

We consider a model where the degenerate 3c? electrons 
of the active site are coupled to the neighbouring oxygen 
orbitals by a hopping term modulated by Slater-Koster 
parameters. The oxygen orbitals are on their turn cou- 
pled to the 3d orbitals of the inactive Mn*+ sites. The 
Hamiltonian of the total system is 

H = Ha+Ti+T2 + Hi+H2, (1) 

where Ha is the atomic Hamiltonian for the active Mn 
site and Hi and H2 are the Hamiltonians for the first- 
neighbour oxygen atoms and the neighbouring inactive 
Mn^+ sites, respectively, in the absence of hopping. The 
Ti and T2 are the hopping terms for the Mn'^+-0 and 
0-Mn"'"'" bonds, respectively. 

Ti = V2tY^ {gol,^d.^2 ,^ + ol,^d^2 ,^ + ol,^dy2 ,^) +c.c., 

(2) 



where the z and x axes lie in the Mn02 plane, t is the 
Slater-Koster parameter, g is the reduction factor of 
the hopping along x (which is taken parallel to the zig- 
zag segment). The d and o are the second-quantization 
operators for Mn 3d and O 2p electrons, respectively. Ne- 
glecting the smaller 14 parameter, we consider only three 
oxygen orbitals (six including spin) and o^.y^z represents 
the orbitals in the x,y,z directions. Only one oxygen per 
direction is considered and a factor \/2 is included in 
Ti, so that the o's represent symmetrized orbitals. The 
, dy2 and dz2 are linear combinations of eg operators 
pointing along the three cartesian directions (e.g., dx2 = 
y/3/2 dx2-y2 + l/2dz2). The Oz and Ox degrees of freedom 
are on their turn coupled to two Mn"*"*" sites by a hopping 
term 

T2 = tY^ {ol.X, + ol^Z,) + c.c, (3) 

a 

where X (Z) represents an eg orbital at the Mn^+ site 
along the X (Z) direction. 

The Hamiltonian Hi for the isolated oxygen atoms is 

^1 =II[epII4^o.,, + t/pp(l-ot^o.T)(l-o,Va)], (4) 

i a 

with i G {x, y, z}. The Hamiltonian H2 is 

H2=Y.[ea + h{]^ + ai)]XlX, + Y,[ed + h{^--ai)]ZlZ,, 

(T a 

(5) 

where h is the exchange energy term that takes into ac- 
count the opposite magnetisation of the two Mn'*+ sites. 
The Hubbard correlation term is absent in Hi^ but in 
the calculation we limit the Hilbert space by disregard- 
ing states with doubly occupied X or Z orbitals. 



3 Choice of parameters 

The parameters for the atomic Hamiltonian Ha-, which are 
the spin-orbit interactions and Slater integrals - except for 
the monopole terms which are strongly screened - have 
been obtained using Cowan's Hartree-Fock code ^31 aver- 
aging the values for the Mn^+ and Mn^+ configurations. 
The atomic parameters are listed in Table 1. We apply 
rescaling factors of 0.75 and 0.8 to the dd and ■pd Slater 
integrals, respectively, as is common practise to obtain ef- 
fective Slater integrals for correlated transition metal sys- 
tems 11 . For we start from the estimated Udd of 
about 5 eV from previous studies based on spectro- 
scopical data. Considering a starting configuration i2gj,^gj. 
and adding two extra electrons we obtain i2gX^si^2gT " There- 
fore the experimental Udd is the interaction between a t2g\ 
and an e^x electron, 

Udd = ~ 4/49 F}^ -2/147 F^^, (6) 

from which F^^ is deduced. 

The hopping t is estimated as 1.8 eV and C/pp as 5 
eV. The bare energy Cp of the oxygen orbitals is deduced 



A. Mirone et al.: Spectroscopy of orbital ordering in Lao.5Sri,5Mn04 



3 




640 650 660 

Energy(eV) 



Fig. 2. Experimental diffraction spectrum (dashes with dia- 
monds) and calculation (solid line) for an optimal distortion 
and correlation. 

from the charge-transfer energy A = A eV. This experi- 
mental parameter is defined as the energy to transfer an 
oxygen electron onto a bare Mn atom (in the absence of 
hybridization). Therefore, 

e, = -A + naF2, + GF^, - 14/49 Fj, - 14/49 F^, 
-2/5 Gl, - 9/35 Gl, + {2- np)Upp + ,5hyb, (7) 

where and Up are the average occupation numbers of 
the Mn 3d and O 2p shell, respectively, in the ground state 
of the model Hamiltonian and Jhyb is the residual energy 
shift of the oxygen orbitals when we consider hybridization 
with inactive sites only. The value that we find is ep — 54.9 
eV. Of course, such a large energy value is not referenced 
to the vacuum level but to a bare 3d orbital with zero 
occupation. The spin-magnetisation parameter h for the 
Mn'*+ sites is 2.5 eV in the case of full spin ordering. The 
energy of the bare orbitals X and Z is 

ed^ep + A + Sd, (8) 

where Sd is a free parameter that is allowed to take val- 
ues in the order of 1 eV, which accounts for different ef- 
fects, such as the charge splitting between the sites and 
an additional effective ligand field that compensates for 
the incomplete inclusion of hybridization for the X and 
Z orbitals in the model. The spin quantization axis z in 
Eq. which fixes the orientation of the magnetisation for 
the Mn'*+ sites, lies in the (001) plane and we have tested 
several different orientations. 



4 Results 

Figure 2 shows the calculated spectrum for a hybridiza- 
tion reduction factor g = 0.7 and a Mn*+ energy term 
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Fig. 3. Calculated diffraction spectra for different values of 
the distortion and magnetic correlation parameter. Solid line 
is for g = 0.7; h = 2.5. Dashed line is for g — 0.85; h — 1.25. 



5d = 0.8 eV. In the calculation the Hilbert space of the 
ground-state configuration and the excited-state configu- 
ration (with a 2p hole and an extra electron in the valence 
band) have been fully expanded with the constraint to 
have between 8 and 12 electrons on the Mn^+ sites and 
zero or one electron in each of the X, Z orbitals. The 
scattering factors are calculated in the dipole approxima- 
tion and using a Lorentzian line width broadening of 0.45 
(0.55) eV for the {L2) edge. The magnetisation is taken 
along the [110] direction which gives the best fit. The re- 
flectivity is calculated as the squared 00 scattering fac- 
tor divided by the absorption of the sample [p. The ab- 
sorption is calculated as a stoichiometric average of the 
Lao.sSri 5Mn04 elemental absorption. The absorption for 
Mn is obtained using the experimental absorption mea- 
sured for Lao.5Sri.5Mn04 P at the ^2,3 absorption edges 
and joining it with the tabulated values. The experimental 
spectrum has been measured at 134 K pQ. 

Our calculation reproduces all the spectral features, 
the only problem is that the right-hand shoulder of the 
lowest energy peak has too much intensity and appears as 
a separate peak. It is interesting to look at the behavior 
of the system as a function of the oxygen displacement. 
Going from zero distortion (g ~ 1) to strong distortion 
{g < 0.5) we have observed that the L3 peak is already 
dominant at very low distortion {g ~ 1), while at lower 
g values it loses intensity compared to the L2 peak. This 
result is very different from all previous analyses based 
on simple ligand-field models |14| . which predict a direct 
relation between the L3 main peak intensity and the Jahn- 
Teller distortion. We think that the ligand-field model may 
fail in this case because the resonant diffraction is sensi- 
tive to the scattering factor differences for different po- 
larisations and these differences depend not only on one- 
particle energy shifts but also ( thinking in terms of one- 
particle Green functions ) on the spectral weight transfer 
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Configuration 


Pdd 


Fdd 


Cm 


Fpd 


Gpd 


'-'pd 


Cap 


^dd 


^pd 


Ground state 


8.94 


5.62 


0.05f 


5.86 


4.38 


2.5 




5.7 


^■^F'^d 


Excited state 


9.53 


5.98 


0.063 


5.86 


4.38 


2.5 


6.85 


5.7 





Table 1. Calculated atomic Ifartree-Fock values |13| for the configuration averaged Slater integrals and spin-orbit interactions 
(in eV). 



due to hybridisation, which ligand-field model neglects. 
We think that, as a consequence, ligand-field model cal- 
culations tend to give stronger spUtting of id orbitals to 
compensate for the missing spectral weight effect, in par- 
ticular Stojic et. al f find a 1.35eV splitting for the t2g 
band. Such a value seems unusual to us. 

We obtained the optimal fit for g = 0.7. For g val- 
ues between 1 and 0.6 the magnetisation of the Mn^+ 
site is parallel to that of the Mn^+ site which lies along 
the X direction. For g values below 0.6 the magnetisa- 
tion is reversed. What happens is that the Cg orbital, 
laying along the 3x^—r^, ahgns its spin parallel to that 
of the Mn^+ atom which is found in the Z direction be- 
cause for too small g the effective hopping on the X or- 
bital becomes weaker than on the Z orbital. However this 
phenomenon could be an artefact of our model. In fact, 
the model neglects the t2g hybridization which always 
gives an antiferromagnetic coupling with the Z orbital. In 
any case beyond this, possibly artificial, magnetisation- 
reversal threshold the L3 main peak is strongly enhanced 
compared to the L2 and it shifts to lower energy by about 
1.5 eV. This behavior resembles what happens in the case 
of the Pro.6Cao.4Mn03 system below Tn [IS . It is interest- 
ing to test with the model to see what happens when the 
spin correlation is reduced. Our model, in particular the 
H2 term in Eq. considers a well established spin order 
and is simplified by using symmetrized orbitals. A way to 
mimic the reduction of the spin correlation, without losing 
the simplicity of the model, is reducing the value of h in 
Eq.O We show in Figure 3, with the solid line, the diffrac- 
tion spectra for reduced values of the distortion {g = 0.85) 
and spin correlation {h = 1.25 eV). The solid line is the 
optimal fit to the experimental data (from Figure 2). The 
change in peak heights reproduces the experimental be- 
havior between Tq and Tn PP, i.e., the main L3 peak loses 
intensity and the right and left shoulders of the L2 gain 
in intensity. 

Analysing the one-particle Green function for the ground 
state we find that in our model the Cg electron spends 
only 15% of its time on the Mn''+ site. This is a very pro- 
nounced charge separation. We also find that the occupied 
eg orbital at the Mn^+ site has 3x^—r^ symmetry. This 
result contradicts the interpretation of the experimental 
x-ray linear dichroism at the Mn ^2,3 edges by Huang 
et al. |16| . They have observed a stronger absorption for 
in-plane linear polarisation than for out-of-plane polarisa- 
tion, from which they deduce that the occupied Cg orbital 
at the active sites has x'^ — symmetry. 



5 Conclusion 

We have analysed the experimental soft x-ray resonant 
diffraction of the orbital ordering in the half-doped man- 
ganite Lao.5Sri.5Mn04 using a many-body cluster calcu- 
lation. The dependence on the order parameters of the 
peak intensities calculated by our model is very differ- 
ent - if not completely opposite - to those found using 
simpler ligand-field models. In our model a central Mn^+ 
site hybridizes with the first shell of neighbouring oxygen 
atoms. On their turn the oxygen atoms are hybridized 
with neighbouring Mn''+ atoms in the Mn02 plane. The 
experimental spectrum at 134 K has been reproduced with 
good agreement, using a slight distortion of the planar 
oxygens and by supposing a strong local magnetic cor- 
relation between Mn sites. The temperature dependence 
of the spectra between Tq and Tn has been reproduced 
by reducing simultaneously the oxygen distortion and the 
spin correlation. It is not possible to reproduce the tem- 
perature dependence of the spectra between Tc and Tn 
by simply reducing the distortion. This shows that spin 
correlation is an important ingredient between Tc and Tn 
and increases with the lattice distortion. We also found a 
pronounced charge separation between the two sites and 
an in-plane orbital polarisation. The latter does not yet 
allow us to explain the the results of some recent x-ray 
linear dichroism observations. 

We thank the ESRF computing service for computa- 
tion of the spectra, in particular Wolf Dieter Klotz and 
Gaby Forstner for taking care of computer clusters and 
grid computing. 
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